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In the multi-fields Phantom power law we investigate the analytical behavior of many 
scalar helds working collectively, where ipi is the ith scalar field. Furthermore, we evaluate its 
parameter values by applying certain constraints on our model parameters, and then compare 
these values with current observational data from Observational data of CMB, BAO and Hq. 
1— ij Through our results, we deeply observe that in the dark-energy EOS parameter at the Big 

Rip always finite with the pressure and dark energy density divergence. 
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CM ■ 1 Introduction 

>: 

Cosmological observations made, by end of the last century and beginning of this century, have 
conclusive evidence for late cosmic acceleration [1J. It is driven by an unknown fluid violating 
strong energy condition (SEC), such that oj < — ^, with oj being the ratio of pressure density 
to and energy density. This exotic fluid is known as dark energy differential equation (DE). 
Recently, Caldwell argued that experimental observations favor the case oj < — 1 more than the 
case oj > — 1 also disobeying the weak energy condition. Also the Dark energy, obeying oj < — 1, 
is dubbed as phantom [2]. With the use of General Relativity (GR)which based on Friedmann 
equations, it observed that phantom-dominated epoch of the universe goes faster, but ends up 
in the form of big-rip singularity in a finite future time [3]. Phantom dark energy fields are 
characterized by violating the main energy condition, p + p > 0. Also the conservation equation 
has the striking consequence that the energy density increases with expansion. When p = up 
, with oj < —1 , are those condition when matter called Phantom energy [3J. On the basis of 
observational data R.R. Caldwell noted that EOS parameter "oj" has a very, short range near 
oj = — 1 with more likelihood to the side of u; < —1. He argued that this possibility could not 
be neglected for the dark energy fluid. The scale satisfying the condition oj < — 1 is known 
as phantom. Alternatively, a very good description of the evolution of universe, is discussed in 
Reference [5j [6] and [7] . 

As energy density and pressure for a homogeneous scalar field is given by p^ = \(p 2 + V((p) 
and p^ = ^(p 2 — V(<f). For phantom scalar field ip , w = using the condition w < —1 in 
above equation, a surprising result is obtained that (p 2 < 0, showing negative kinetic energy for 
the phantom scalar field. A power law cosmology is defined by the cosmological scale factor 
evolving as [8]. 
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We want to bring down observational constraints on single field to extend this rule for Nfields 
power law. It is very interesting to generalize the above studies to multi-fields phantom power law 
with other potentials. For the simplicity of our model we take the isotropic and homogenous 
and Friedmann-Robertson- Walker (FRW) metric having curvature "fc" and "t" is the cosmic 
proper time. In this paper, we study a general Nfiels power law cosmology with the scale factor 
given in terms parameter (3 with out any dimension [8] . In the case of Multi- fields , we modified 
the scale factor a(t) = ao((t s — t)/(t s — to)) 13 in order to achieve the self-stability, where t s is a 
required positive reference time [9]. 

The field starts from near an unstable equilibrium (taken to be at the origin) and climb up 
the potential to a stable maxima. The Phantom model in which the observable Big Rip occur 
during the climb up of scalar field and its magnitude is at most of order M p . 

The arrangement of remaining sections of this paper is as; Section 2, we formulate the whole 
picture of phantom power-law cosmology for multi-fields. Section 3 we consider observational 
data to impose certain bounds on the Multi-fields parameters, and finally, Section 4 is devoted 
for summery and discussion. 

2 Nfields with Power-Law Expansion 

In this section we present Phantom cosmology under power law expansion, when many fields 
are working collectively with <pi is the ith phantom scalar field. We assume the homogenous and 
isotropic Friedmann-Robertson Walker (FRW) background metric, 



ds 2 



dt 2 - a 2 (t) 



dr 2 
1 — kr 2 



+ r 2 dn 2 2 



(1) 



Where is 2-dimension unit sphere volume, t is the cosmic time and k represents the curvature 
of 3-dimensional space with k = 1,-1,0 corresponds to open, flat and closed universe respec- 
tively. 

Our model is given by the following action, reads ([10j.[TT]): 



s = j d^xv^Yl 



(2) 



Where R the Ricci scalar, Vi(4>i) is potential of phantom field and G the Newton gravitational 
constant. The Lagrangian L m stands for the total matter of the universe including (dark plus 
baryonic). Finally, we concentrate on small redshifts, therefore, we are neglected the radiation 
sector, with speed of light as unity [12]. 



Through out our work we assume the flat geometry of the universe i.e. k = 0, for multi-field, 
the Friedmann equations write, reads ([1~3].|14|) 



H2 ~ ^E^ + ^l' 



(3) 



where H(t) = a/a is the Hubble parameter represents the expansion rate of the universe at time 
t and M p = (8ttG)~2 is the Plank mass. In the above expression p ( p i and are respectively the 
energy density and pressure of the ith phantom field which is denoted by 4>i , which are given by 
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h = = -^£<^ 2 + £ wo- 

i i i 



(4) 



P<t> 



(5) 



Moreover, the evaluation equation for the multi-field is defined as 



£ Pfa + 3fT £(p^ + - 0. 

i i 

After simplification, the above equation becomes, 



(6) 



£^ + 3ff£<^-£-^o. 



(7) 



As in phantom cosmology the dark energy sector is attributed to the phantom fields, and thus 
its equation-of-state parameter is given by 



UJ, 



Prrij 
Si Pm,i 



(8) 



where w m = ]T\ ^m, • 
For the matter density 



£p mi + 3#£(1 +Vm i )Pm i ^ 0, 



with the simple form of its solution is 

£^ = £^f 



(9) 



(10) 



where n« = 3(1 +u mi ) and we are discussing only massive scalar fields, but the case of massless 
scalar fields are neglected here. As a special case when uj mi is equal to zero, then n = 3. 
Also from Eq. (Hj) 



H 



1 



6M| 



3 £^-£ 



Prrii fH 



(11) 



When we are working with Phantom cosmology, we replace t by t s — t; the reference time t s 
is sufficiently positive, then we obtain 



a(t) = a 



t s -t 







,t s — to, 

with the Hubble parameter and its derivatives with respect to time is 



(12) 



3 



m K (13) 

Thus for the value of (3 less than zero, we find an accelerating "a(i)" is greater than zero 
and expanding "d(i)" is greater than zero, universe, we find that H (t) is positive therein, which 
implies that it provides super acceleration, this is only possible for phantom power-law cosmol- 
ogy. In addition, for the exponent f3 < 0, and at late time t = t s , the scale factor a(t) and H(t) 
of the universe both are diverge, as a result it goes to a Big Rip. 

But then, the important point which is already discussed in reference |15| . We find the 
potential J2i V i(t), by using Eqs.© and (0} in Eq. (HI]). Since n; = 3(1 +w m j; 1 < i < N. 
For dust dominated universe uj m i —> 0, this implies that re, = n = 3 and as a result we obtain 
Y2 n i = 3N. Therefore 



V(t) = J2Vi(t) = N 

i 

From Eq. (jlip we can obtain 

4> 2 (t) = N 



Mi 



3f3 2 -p 



hp 



mo 



6 a J 



t s — to 
U-t 



+ 



t,-t 



where 0(t) =Yl l %^i- 

Using the values from Eqs. (fT4"|) and (fl~5j) in Eq. ([!]) 



3/3' 



(14) 



(15) 



3/3 2 4p mo f t s -t 
(t s -t) 2 3a 3 V^-t 

Again putting Eqs. (fl~4"|) and (fT5j) into Eq.([5|) we get 



3/3 



(16) 



P0 ~ TV 



^(-3/^-/3 + 3/3) _ _ + Pm0 ft. -to 



(t - t s f K SJ w 3a 3 V *s - * 

In the case of phantom Nfields, the dark energy equation-of-state is 



UDE{t) =V4>IP4> 



which implies that 



3/3' 



(17) 



w DE {t) = {-! + -) 
We see that for Big Rip behavior uj£,E{t) always having finite value [T6] , 



(18) 



For /3 less than zero possesses additionally a positive H(t) that leads to super-acceleration 
|17] . So such kind of scenario expansion is always came with acceleration. 

Furthermore, with the value of j3 less than zero, at t equal to t s the scale factor and the Hubble 
parameter diverges, that is the universe results to a Big Rip. We investigate that these behaviors 
of Nfields phantom cosmology are very similar with result of single field phantom cosmology with 
power-law [15]. 
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3 Observational Bound For Nfields With Model Parameters 



In this section, we apply the techniques that conform the values of multi-fields with the observa- 
tional data. Now we are fitting the observational data, presenting our results in the case of many 
fields work collectively. We observed that our values for all parameters in Table I are best fitted 
with minor error of accuracy, which is negligible for large scale, we also provide the la bound 
of every parameter. Similarly in Table II we present the maximum possibility of the values up 
to la bound for the derived parameters, namely the power-law exponent /3, the present matter 
energy density value p mo , the present critical energy density value p CQ and the Big Rip time t s . 
As we observed that f3 is always less than zero, as expected in consistent phantom cosmology. 
Furthermore, we observed that the BR time is one order of magnitude greater than the present 
age of the universe, which shows that such an outcome is predictable in phantom cosmology, 
unless one include additional mechanisms as shown in 1 191. 



V(t) « N 



6.5 x 10 27 



2.52 x lQ- 67i x (3.30 x 10 i8 - t) 



vl9.54 



(3.30 x 10 18 -t) 2 
While when we consider WMAP7 data alone, it provides 



(19) 



V(t) « N 



6.41 x 10 



27 



(3.30 x 10 18 - t) 



1.98 x 10 _db9 x (3.30 x 10 i8 - t) 



vl9.37 



(20) 



Here we noted that in above results although the second term is very small at early times of 
the universe, but becomes very important at late times, this situation is close to the BR. Now 



the scalar field evolution at late time (t 
Then we obtain the new result 



t s ), Pm can be neglected and also set ao = 1. 



<p(t) 




(ts-t) 2 



dt. 



(21) 



Furthermore, for the combined WMAP7 + BAO + Hq, we get 



<p(t)^-^j2Nf3M$ln\t s -t\, 

which implies that 

<p(t) ^ \/iV[-2.645 x 10 13 * In |3.30 x 10 18 - t\ 
while single WMAP7 data provides that 



(22) 



(23) 



<p[t) = \/iV[(-2.643 x 10 13 ) * In |3.30 x 10 18 - t\] (24) 



According to our exactions the phantom field and the kinetic energy diverges at the Big Rip. 



Parameter 


H + WMAP7 + BAO 


WMAP7 


to 

#0 

^CDMO 


13.78 ± 0.11Gyr[(4.33 ± 0.04) x 10 17 sec] 
70.2~}j km/s/Mpc 
0.0455 ± 0.0016 
0.227 ±0.014 


13.71 ± 0.13Gyr[(4.32 ± 0.04) x 10 17 sec] 
71.4 ±2.5 km/s/Mpc 
0.0445 ± 0.0028 
0.217 ±0.026 
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TABLE I: Observational maximum likelihood values in la confidence level |18| . 



Parameter 


H + WMAP7 + BAO 


WMAP7 


PmO 
PcO 
ts 


-6 ^2- u - 2b 

°- d ^+0.24 

(-6.532 ± 0.38) x 10~ 27 kg/m 3 
{9A-° A 3 ) x 10~ 27 kg/m 3 
104.83+J;? Gyr [(3.30 ± 0.06) x 10 18 sec] 


-6.51 ±0.4 
(2.513 ± 0.27) x 10" 27 kg/m 3 
(9.63 ±0.58) x 10~ 27 kg/m 3 
102.8 ± 3.4 Gyr [(3.24 ± 0.21) x 10 18 sec] 



TABLE II: Derived maximum likelihood values in la confidence level for the power-law 
exponent. 



4 Summary and Discussion 

From Eq. (f23"|) and Eq. (|24p . we find the different values of t and put them in Eq. (|19p anci Eq. 
(|20p . to evaluate the potential for the Nfields Phantom power-law i.e. V (<£>)• 

Thus, for the WMAP7 + BAO + H and WMAP7, we obtain different values of V(<p) 

« iV[6.48 X l 27 e v^(0.77xlO-V) _ 2 52 x 10 - 371 e v / lV(-7.44xlO-V) ] 

and _ 

« iV[6.39 X l 27 e vW(0.78xlO- 1 V) _ 1.9 8 x 10 -369 e V^(-7.5xlO- 1 V) ] 

respectively. 

Through out our work we study the phantom cosmology for multi-fields, we observed that the 
cosmic scale factor a(t) is obeying the power law. When we construct the whole scenario, we fit 
the observationally data of WMAP7+BAO + Hq and WMAP7 alone by applying bound on the 
multi-fields by focusing on exponent (3 and the BR time t s . By using separately WMAP7 data, 
we obtained the value f3 = -6.523 ± 0.38, while the BR is observed at t s = 102.8 ± 3.475 Gyr. 
However, the dark energy equation of state parameter ujde lies below the phantom divide, it was 
expected and at BR time it always remains finite and equal to —1.1533. Although the phantom 
dark anergy density and pressure are diverge at the BR. By using WMAP7 + BAO + H$ data 
set alone we find j3 — 6.51T_g'||, while the Big Rip is observed at t s w 104.5^ 2 'g Gyr, in la 
confidence level. Definitely, the subject of Nfields quantization of such scenarios is open and 
needs further investigation on this. 
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